The Lasting Effect of Initial Conditions on Single File Diffusion 
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We study the dynamics of a tagged particle in an environment of point Brownian particles with 
hard-core interactions in an infinite one dimensional channel (a single-file model). In particular we 
examine the influence of initial conditions on the dynamic of the tagged particle. We compare two 
initial conditions: equal distances between particles and uniform density distribution. The effect is 
shown by the differences of mean-square-displacement and correlation function for the two ensembles 
of initial conditions. We discuss the violation of Einstein relation, and its dependence on the initial 
condition, and the difference between time and ensemble averaging. More specifically, using the 
Jepsen line, we will discuss how transport coefficients, like diffusivity, depend on the initial state. 
Our work shows that initial conditions determine the long time limit of the dynamic, and in this 
sense the system never forgets its initial state in complete contrast with thermal systems (i.e a closed 
system which attains equilibrium independent of the initial state). 

PACS numbers: 05.40.Jc,02.50.-r 
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I. INTRODUCTION 

Particles diffusion in one dimensional systems, with 
hard core interactions, have been studied for many years 
|l|-[l6|. One aspect of this problem is the motion of a 
tagged particle. This kind of system can be used as a 
model for the motion of a single molecule in a crowed 
one dimensional environment such as a biological pore 
or channel [13, [l3, and experimental studies of physi- 
cal systems such as zeolites [13] and collo id p articles in 
confined topology [20] or optical tweezers |2]| . 

In a system of interacting Brownian particles the mo- 
tion of a tagged particle was thoroughly investigated. 
The initial state of the system of particles is in most 
previous works taken from equilibrium. That means par- 
ticles are initially uniformly distributed with density p. 
One of the known results is that the tracer particle subd- 
iffuses, i.e. (a;f-) ~ 2I?i/2i^'^- Harris was the first to pro- 
vide a theoretical derivation for this phenomena by using 
statistical arguments [l|. For a finite system, e.g single- 
file diffusion in a box, the tagged particle's mean- square- 
displacement reaches equilibrium, limt_j.oo i^'r) ~ const. 
[23 . l23j (see also [2J] for periodic boundary condition) . 

An intriguing treasure is found in an Appendix of 
Lizana et al. [25j. They note by passing that the gen- 
eralized diffusion coefficient -D1/2 is sensitive to the way 
the system is prepared. That is a surprising result since 
we expect diffusivities of interacting system not to be 
sensitive to the initial conditions. 

Here we confirm this prediction using the Jepsen line, 
showing that the prediction in [2^ is correct. To show 
that Di/2 is sensitive to the initial state of the system we 
consider two ensembles of initial conditions: an initial 
state that is taken from equilibrium versus particles ini- 
tially situated in equal distances between each other. We 
show that Di/2 for both ensembles differs by a prefactor 
of V2. 

Since transport coefficient like l?i/2, for the many- 
body interacting system, is sensitive to initial prepara- 



tion we must challenge basic concepts in non-equilibrium 
statistical mechanics. For example in sec. V in investi- 
gate the response of the system to external field. Will 
that response depend on initial condition? and how do 
we formulate the Einstein relation, if at least diffusivity 
Di/2 is sensitive to the initial condition. Further we in- 
vestigate the correlation function {x{t + A)x{t)) of the 
process showing clear differences between the ensembles 
of initial conditions (see sec. HI- IV). These are found also 
in long time limit A — >■ cx) where the correlation function 
is markably different from the mean-square-displacement 
(MSD). 

Finally we investigate the time average MSD (sec. VI). 
Will that time average depend on the initial condition 
(like the ensemble average) becomes an interesting ques- 
tion. Further is the MSD ergodic in the sense that for 
two identical initial condition do the corresponding time 
and ensemble average MSD coincide? 



II. THE MODEL 

In our model we have 2N + 1 point identical Brownian 
particles, with hard core interactions, in a one dimen- 
sional system. D represents the diffusion coefficient of a 
free Brownian particle. We tag the central particle, so 
there is an equal number of particles on its right and on 
its left. Initially the tagged particle is situated at the 
origin xxit = 0) = ■ The system is stretched from —L 
to L. The size of the system and the number of particles 
is infinite {N,L — > (X)), but N/L = p is fixed, where 
p represents the particles density. p~^ is the mean dis- 
tance between nearest neighbors. We consider two types 
of initial conditions. The first is the case of particles 
distributed with fixed density, namely the distance be- 
tween particles are exponentially distributed with mean 
distance between particles is p~^ . This case was treated 
previously. The second case we consider is the case where 
particles are initially situated on a lattice, with equal dis- 
tances between particles L/N = a where a is the lattice 
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FIG. 1: An illustration of the model of particles in a narrow 
channel so particles cannot pass each other with a marked 
tagged particle. The system can be mapped into a chain of 
beads which are interconnected with springs (see the Harmo- 
nization method below). This schematic diagram presents the 
initial state of the system with equally spaced particles. 



constant. We labeled the tagged particle as n = and 
the particles to its right are labeled n = 1,2... accord- 
ing to their order. Similarly, the particles to the left of 
the tagged particle are labeled n = — 1, — 2 . . .. Then the 
initial position of each particle is represented by: 



Xn{t)\t=o = na 



(1) 



where n is the label of the particle, n e {—N... N} (see 
Fig. [J). 



III. THE JEPSEN LINE 

A. Mean-Square-Displacement 

In this section we find the mean-square-displacement 
(MSD) with a rigorous method by treating the problem 
with a theoretical tool called the Jepsen line (see details 
below) 0, i, H, m. We find the MSD of the tagged 
particle: 



{xU't))iat =aV2J—Vt 



{xUt))um='2p-\j-Vi, 
TT 



(2) 



which is in agreement with the statement in ^25]. Here 
{.■.)iat refers to ensemble average when the system is ini- 
tially on a lattice (non-equilibrium case), and {.■.)uni 
refers to an equiUbrium state (uniform distribution). 
Notice that the difference between the mean-square- 
displacement (MSD) of the two initial conditions is a pre- 
factor of %/2 (see figE]). More specifically we compared 
between the two initial condition by taking a = p^^ since 
that is the average separation between particles on the 
lattice. This is a surprising effect, since we find that the 
influence of the initial conditions on the diffusion of the 
tagged particle in infinite system is lasting forever. 
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FIG. 2: The MSD of the tagged particle in two cases: When 
the particles are initially in thermodynamic equilibrium (up- 
per line), and when the particles are initially with an equal 
distance between each other (see simulations details in ap- 
pendix A). Comparing the two results gives the factor of \/2 
(see Eq. ([2])). The averaged distance between particles is 
a — p~^ . 



B. Jepsen Line Method 

The motion of a single particle without interactions 
with other particles is given by the Green function 
g{x,Xo,t), where g(x,xo,t)dx is the probability of the 
non-interacting particle which started at Xo to be in 
(x, X + dx) at time t . For an infinite system the Green 
function for a free Brownian particle is simply a Gaussian 



l{x,xo]t) 



1 



^/A■KDt 



(3) 



We consider initial positions of particles given by 
Xn{t)\t=o = an where ne{—N. . . N} is the label of the 
interacting particles and a is the lattice constant. A 
straight line that is initially at a; = and follows x — vt 
{v is a "test" velocity) is called the Jepsen line (see fig. [3]). 
Initially there are iV -I- 1 particles, including the tagged 
particle, to the right of the line, and N particles to the 
left of the line. 

Let a be the label of the first particle that is situated to 
the right of the Jepsen line, therefore initially a(t)\t=o — 
0. The random variable a increases or decreases by -|-1 or 
— 1, according to the following rules: if a particle crosses 
the Jepsen line from left to right a decreases: a ^ a — 
1 and if a particle crosses the line from right to left a 
increases: a — > a+1. Thus a is a random walk decreasing 
or increasing its value -1-1 or —1 at random time (see fig. 

n. 

In one dimension, a hard core elastic collision event, 
the result of two identical particles (same mass) colliding 
is that they switch their velocities. For over damped 
Brownian particle this is equivalent to two particles that 
pass through each other, and after the particles cross 
each other, the labels of the two particles are switched. 
Instead of relabeling the particles after every collision, we 
let particles pass through each other, and then at time t. 




We define 



FIG. 3: A schematic example of the trajectories of particles. 
The straight solid line that starts at (0, 0) and follows x = vt 
is called the Jepsen line. We label the particles so that the 
tagged particle is n = and the particles to its right are 
labeled n = 1,2... according to their order while the particles 
to its left are labeled similarly with n — —1, —2, . . . 



we label our particles. So in fact in the interval (0,i) we 
view the particles as non-interacting (see fig. [S|) 



N 



E^^ 



(4) 



where (5a„ — a^ 



>;f^„ and a^ is the number of times 
that nth particle crossed the Jepsen line from right to left, 
minus the number of times that it crossed the line from 
left to right and a^„ is defined similarly. For example 
a{t) = means that all the particles stay on their original 
side of the Jepsen line, so the first particle that is situated 
to the right of the line, is the tagged one. The variable 
a is determined by 2N + 1 random variables and since 
iV >> 1 we can neglect the contribution from SaQ. 

For calculating the probability density function PAr(a), 
we mark P^ (a^o") as the probability that the n-th particle 
starts at a;o to the i side of the Jepsen line and ends at 
the j side of the line {i,j G {R,L}) [R represents the 
right side of the Jepsen line, and L stands for left] . On 
a lattice xq" = an where a is an equal distance between 
particles. For one step of the random walk, Sa„ can get 
the values —1,0, 1 with the probabilities: 



P((5a„ = 1) = PBL{an)PLL{-an 

5an = ■( P{5an = 0) = PRL{an)PLR{-an 

P{5a.„ = -!) = PRR{an)PLR{~ 



- PRR{an)PLL{-a7i) 



(5) 
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FIG. 4: When a particle crosses the Jepsen line from right 
to left (or left to right a increases) a — >■ a + 1 or decreases, 
a — > Q — 1. 



Notice that 5an depends on the motion of two particles 
initially at an and at —an. For example, if one particle 
starting on an , i.e. right side of the Jepsen line (i?), 
switches to the left {L) of the line while corresponding 
particle on —an (L) remains in L we have 5an = 1. The 
probability Pij [an) is given by Green function and initial 
condition, for example: 



PRL{an) = 



PRR{.an) 



1 (x-a„)^ , 

=e 4Dt dx, 



L Va^kdI 



vt VinDt 



1 _(x-a,i)2 

e ^^^dx. (6) 
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FIG. 5: Illustration of collision event: since the system in 
unidimensional and the particles are identical, when particles 
collide the velocities are switched (B.). Instead of switching 
their velocities, we can switch their labels, and treat each 
particle as a free particle without interaction (A.). The inter- 
actions come into account in sorting the particles and finding 
the central one. 



C. The Motion of the Tagged Particle 

When N -^ oo, a is normally distributed, according to 
the central limit theorem for the random variable a (see 
appendix A), hence 



PNia) 



1 



_(o-(Q>yf 



\/27rcr2 



(7) 



where (a) is the mean of a and a'^ = (cP'\ — (a) is its 
variance. 

The probabihty to find the tagged particle to the left 
of the Jepsen line is: 



/vt 
Pt{xT)dxT, 
-OQ 



(8) 



and our goal is to find the PDF Pt{xT) of finding the 
tagged particle in {xt,xt + dxr)- The event xt < vt 
is equivalent to the case that the first particle that is 
situated to the right of the Jepsen line is labeled a > 0, 
hence 



P{xt < vt) 



1 



(a-(°>)^ 



\/2tt(j^ 



'da. 



(9) 



Taking the derivative of ©, following Eq. ([5]) replacing 
vt ^^ XT, we find the probability density function (PDF) 
of the tagged particle position 



Pt,N{xT) 



V2^ 



1 _M^ 

= 6 2<T^ 



(10) 



where we will soon take the large N limit. By definition 
the average of variable (5a„ is : 

(fc„)= Y. Sa,-,P{San). (11) 

By using PLLi~an) = 1 - PLR{-an) and PuR^an) = 
1 — PjiL{an) it is easy to see 



((Sa„) = PRL{an) - PLR{-an). 



(12) 



Notice that here Prr, Plr etc. (Eq. ^) is now calcu- 
lated with XT is the upper or lower integration bound, 
i.e. replacing vt — >■ xt in Eq. ([6|). Similarly we can find: 

{Sal) ^ PRR{an)PLR{-an) + PLL{~an)PRL{an). (13) 

For the average we clearly have 



N 



y^ (5a„ 



(14) 



For the variance 



TV N N 

{a') = Y. ('5"«) + E E (^«n^«™) ■ (15) 

n=l n=l m^n 

Note the random variables {(5q„} are indepen- 
dent, therefore the covariance term vanishes, i.e. 
Z]n=i J2m,^n (SanSam) = 0. Since iV ^- oo we can ap- 
proximate Pt.N{xT) near its saddle point. The equation 
for finding the saddle point Xg is: 



E^«^(^ 



PLR{~an) 







(16) 




FIG. 6: The tagged particle location with lattice initial con- 
dition is distributed normally when the system contains large 
numbers of particles. The solid curve is the normal distribu- 
tion with mean zero and variance the Eq. (|27|l . We simulated 
two different times: f = 10 (o) and t = 100 (*). 



the particles are unbiased, therefore Xg — 0. Approxi- 
mating (a) near the saddle point and using Eqs. ^ and 
p2| give (the first non-zero term): 

^ 2 („„)2 

(a) « d^ (a) \ XT = V /-— - e ^°' xt- (17) 

Similar approximation for a^ by using Eqs. ([6|) and (|13|) 
gives: 



2 ^, / 2\ 



(18) 



N 



2E 



■;j V47rL>t 
^ 



1 (:.-r.o)" , 

le 'lo* da; 



1 _xx+ii^ , 
e 4Dt dx. 



V47rDi 



Therefore xt is normally distributed with mean 

(xt) = (19) 

which is clear from symmetry, and variance 



„_i . \^ e -iDt dx ,/ ^ e lut da; 



2 > 1 , e -iDt 



(20) 



D. Mean-Square-Displacement 

We evaluate the MSD of the tagged particle that is 
given in Eq. (f20| . 



N 



1 



1 



^ 1 _(^^ 

lim > — p 'lo* ^ , 

where we used Euler-Maclouren Formula: 

N~l ,N 



(21) 



E/fc= / /(fc)dfc--[/(0) + /(iV)]. (22) 
fc=i -^0 



For the time when particles interact with each other, i.e 
a^/D «t: 



E 



1 (g^)^ 



^ ViTTDt 

Similarly we find the numerator: 



1 
2^' 



(23) 



^ OO 



an , „ „ , an , 1 . / 2 vADt 1 



4^ V4i5t V4I?t 4 

n— 1 ^ ^ 



Therefore, when particles interact with each other. 



a 

-ry = Tint «t 



2' 
(24) 

(25) 



the numerator in Eq. ([20| is 

1 x"^ „ „ , cm. , „ „ , an 
J2Erfci^==)Erfci- 



1 / 2 V4:Dt 



n=l 

hence we find 



/4Dt 



ADt 



4 V TT a 



{^T')lat = \[l^^D-t. 



(26) 



(27) 



For particles that are initially uniformally distributed it 
is well known that [3J: 



\ / UTl't h^' 



(28) 



Comparing Eq. ((27|) and ([28|) we see that if we assign 
a = p^^ the two results differ by a prefactor \f2 . To 
summarize in a long time limit the PDF for the tagged 
particle interacting with a bath initially on a lattice is 



Pt(xT) 



1 



exp(- 



y/av^SnDt aVSDt 



:)' 



(29) 



see fig. ini 



IV. CORRELATION FUNCTION 

A. Correlation Function 

The correlation function between the location of the 
tagged particle at time t to its location at time t + A 
is defined as: {xxit + A)xT{t)) . We will soon show that 
the correlation function for particles initially on a lattice 



{xT{t + A)xT{t))i^t = a\ —{V2tTA - VA), (30) 

V TT 

while for particles initially in an equilibrium state 



{XT{t + A)xT(t))„„ = p-\/-{Vt + A + Vi - VA). 

(31) 



TT 




FIG. 7: The normalized correlation function shows a sig- 
nificant difference between three cases: free particle (upper 
curve), uniform initial condition (middle curve) and lattice 
initial condition (lower curve). When A » t the normalized 
correlation function in thermal equilibrium goes to 1/2 and 
in the non thermal case it decays to zero. The upper curve 
represents the normalized correlation function for a free par- 
ticle. The circles represent the simulations, and the solid lines 
represent the theory Eq. ((33)). 



Which agrees with Eqs. ([271) and ^ when A = 0. 
Ea. (l5T|) has the structure of the correlation function of 
fractional Brownian motion 28] . 

We define normalized correlation function: 



{xT{t + A)xT{t)) 

(xUt)) 



(32) 



Hence using Eqs. dH]), dH]), ([30]) and dST]) we find 

fA 



^" ( T ) = 



A 
1+2^ 



1 



(1 



A 

T 



(33) 



For free Brownian particle the normalized correlation 



function is: 



9 



free 



free ' 



(xHt)) 

normalized g 
A/t — > OO we get : 



1 since {x{t + A)x{t)) ^^ 



When A/t = the correlation function is 
1 (see fig. El). When 



lat 



9" 



^free 



9 



glat 

.uni 





1 

2 



(34) 



This result emphasizes the strong dependence of initial 
conditions at long time limit. 



B. The Harmonization Method 

Up until now we considered hard core interactions be- 
tween particles, and used the Jepsen line to find the 
MSD. To find the two times correlation function we follow 



Lizana et al. [25| and approximate a general interactions 
into harmonic interactions between particles. 

The Langevin equation that describes the motion of 
the n-th particle is: 



? 



dXn (t) 

dt 



J2 F{Xn{t) - Xn'{t)] + Vn{t) + fn{t) (35) 



where Xn (t) is the position of the n-th particle at time t 
and ^ is the friction constant 



D = 



knT 



(36) 



F 



dV 

dXn 



is the force due to the interaction between 



particles, where the interaction potential between the n- 
th particle and the n'-th particle is V{\xn{t) — Xn'{t)\) 
and the potential V has singularity when x„ (t) — Xn' (t) 
so particles cannot pass each other, rjn is white Gaus- 
sian noise with mean zero and covariance {rjn{t)rjn' {t')) = 
2^kBTS{t — t')Snn', and /„ is an external force. 

The main idea of the method is mapping the many- 
body problem into a solvable harmonic chain, by map- 
ping the system into a system with beads interconnect- 
ing with harmonic springs (see fig. [1|). We convert the 
general interactions term in Eq. (I35p to forces from the 
nearest neighbors spring coupling, i.e. 

fir (f) 

^^^ = KK+i(i) + a;„-i(t) - 2x„(t)] + 7?„(t) + /„(t) 

.(37) 
where the effective spring constant k for one dimensional 
point like particles with hard core interactions is j25| 



N^keT 
L2 ■ 



(38) 



In [25J K is given in terms of the compressibility, which 
that means we can treat a general type of interaction 
(beyond hard core). Eq. p7l) is the Edwards- Wilkinson 
equation, whose relation to single file diffusion was un- 
covered in [23] . Under the assumption that the particles 
interact with others 



Taking the Fourier {n —> q) and Laplace {t — ^ s) trans- 
forms 



x{q,s) 



dn / dte^'9"~^*a;(n,i) (41) 

oo Jo 



of Eq. (go]) gives: 



x(g, .) = ^^(g^^-0) + ^Kg,^)+/fa^). (42) 
E,s -(- Kq^ 
Eq. (l40t is the Rouse chain model which describes poly- 
mer dynamics . In that well known model each monomer 
end is treated as a bead while the interaction between 
the beads is harmonic l30ll. 



Evaluation of the Correlation Function 



We now calculate the correlation function 
(xT(i + A)xT(i)) by calculating {x{q,s)x{q' ,s')) first, 
and then using inverse Fourier and Laplace transforms. 
We assume that there is no external force, f{n,t) — 
and the initial condition is x{n,t = 0) = na, so initially 
particles are on a lattice. 



Using Eq. p2|) gives: 

{x{q, s)x(q\ s')) = Anoise{q, q\ s, s') + Ainit{q, q\ s, s') 

(43) 
where the dependence on the initial condition is: 



, . , ,, ^{x{q,t = Q)x{q',t' = Q)) 
A^n^t {q,q ,s,s) ^ —- ■ -r-- (44) 



{^s + Kq^){^s' + Kq''^) 



and the noise term is: 



t » n. 



DN"^' 



(39) 



we can take the continuum limit and turn a;„(i) into a 
field x{n, t) with the equation: 

^dx(n,t) d'^x(n,t) , ^ ^, , ,,„s 

e L ' ^« ^„; ^ +V{n,t) + f{n,t). (40) 



dt 



d'n? 



^noiseyQi Q i ^i ^ ) 



{vjq, s)v{q' , s')) 
i^s + Kq^){^s' + Kq'^) 



(45) 



Using the inverse Fourier transform J^ ^{ ^^^ ^ } — e "I"' 
and the convolution theorem 



(^ /"OO /"OO I — } — 

{n,n',s,s')^—^ dm dm'e^v ^'"-"le^V V|"'-™'| (^^(to^ t ^ 0)x(to', t = 0)). (46) 

K'y SS J— oo J —oo 



When particles are initially on a lattice 

hat 



{x{m, t = 0)x{m', t = 0));^j = a'^mm' (47) 



(to and m' are integration variables in the spatial do- 



main). We are interested in the tagged particle, so we 
integrate Eq. (P5|) for n = n' = 0. Therefore 



Ainit{n = n' = 0, s, s') = 0, 



(48) 



since f^ e '" "^^m,dm = for n = 0. 

J — oo 

The variance of the Gaussian noise in q and s space is: 



{r]{q,s)r]{q\s')) 
thus we can write: 

Anoise(Q7 Q , S, S ) 



4TTCkBTS{q + q') 



Arr^kBTSiq + q') 



{^s + Kq'^){(^s' + Kg'2)(s + s') ' 
The inverse Laplace transform gives: 



AT:kBT6{q + 9') / _^_^(t+t') ^ ^- -l|t-t'| 



(49) 



(50) 



e 



(51) 



Inverse Fourier transform of Eq. (j5ip gives the covariance 
of the location of the tagged particle, i.e. n = n' = 0, 



{xTit + A)xTit))i^ 



ksT 



(V2i + A 



(52) 



Returning to the original parameters using Eq. psp and 
the Einstein relation Eq. (p6)) gives Eq. ([30]): 



{xT{t + A)xT(i))iat = a\/ — (V2t + A - %/A). (53) 

V TT 

The derivation of the correlation function for the uniform 
case, Eq. ([31]) . is done similarly and we do not include it 
here since it is a straight forward. 



V. VIOLATION OF THE EINSTEIN RELATION 

We would like to examine if the mean displacement of 
the tagged particle (xt(0)_f '^^^'^ ^^^ presence of con- 
stant external force F and the MSD without external 
force <^a;|,(i)) obey the generalized Einstein relation 



{Mt))i 



.(4W) 



F=0 



2kBT 



(54) 



This is a general relation valid within linear response 
theory. However, now that we find sensitivity to ini- 
tial preparation of the system, the generalized Einstein 
relation must be checked 

For constant weak external force /(t) = Fq the mean 
displacement is (see derivation later): 



^lat aFaVDi 



^^^(')^^o - j^fT^ 



(55) 




FIG. 8: Simulations of the mean-displacement of the tagged 
particle when _Fo — 0.5 for the two cases: the uniform dis- 
tribution (*) and the lattice initial conditions (o). The solid 
line represents the theory Eq. H67[) . The mean-displacement 
behaves similarly for the two types of initial conditions as is 
explained in the text. 



and 



° pkuT^jT: 



(56) 



The uniform case was treated previously in [23 • Unlike 
the MSD here we do not see any impact of the initial 
condition on the response {xT{t)). Namely here we get 
the expected result: if we set a = p^^ both results are 
identical (see fig. [5]). 

Comparing the results with Eq. ([54)) gives an ade- 
quacy of Einstein relation for thermal equilibrium initial 
condition (uniform distribution), and violation for the 
non-thermal initial condition (lattice): 



(^T(i))r 



- Fc 



\ lat 



(mwf: = ^0 



2fc^r 

V2kBT 



(57) 



notice the factor \/2. In the derivation of the mean- 
displacement it was assumed that the external force is 
weak, Fq << ^£^^ since we assumed the system is near 
equilibrium. For strong force our theory will fail. 



The Mean-Displacement in Presence of 
External Force 



We now treat the case where a constant force F acts 
on the tagged particle, using the Harmonization method 
of Lizana et al. (2^ . Taking the ensemble average over 
Eq. ([42]) and using the fact that {ri{q, s)) — for white 
Gaussian noise 



{x{q, s)) j: = B,mt {q, s) + Bf {q, s) 



(58) 



where the dependence of the initial state is 



Binit{q,s) 



and the force dependence is 



^s + nq^ 



Bfiq,s) 



^s + nq^ 



(59) 



(60) 



First we evaluate the initial conditions term Binn {q, s). 
The inverse Fourier transformation J^~^ < ^^'f ^ \ = 
g-a|n| g^]-j(j |;]^g convolution theorem give: 



Taking inverse Laplace transform (s -^ t) gives 



{xT[t)) 



(^^^)T7lX —^^^'-^'^^^- (66) 



For constant external force /(f) = Fq the mean displace- 
ment is: 






/tt^k 



(67) 



Btmtin,s) 



- - 

2 V KS 



dme-V »l"~"l(a;(m,i = 0)). 

(61) 
The mean initial position is the same in both initial con- 
ditions, the equilibrium state and the lattice state, i.e 
{x{n,t = 0))„„ = {x{n,t = 0));^j = na (note a = p-^) 



Returning to the original parameters using Eqs. p6p and 
([38| gives: 



{xr{t))t = {-Tit))X - f^- (68) 



dme-V^I^-^lm-a 



V 4ks ^ 

(62) 
for the tagged particle n = therefore the contribution 
from initial state term will vanish: 



B'Zti.^,s) = B^::,{Q,s) 






(63) 
We see that the response is not sensitive to initial condi- 
tions, while the MSD is. 

For the force term Bf{q,s) the inverse Fourier trans- 
form is 



Again here a — p "'^ is the mean spacing between parti- 
cles. 

VI. TIME AVERAGE MSD 



In experiments in many cases we measure an average 
over time [32]. Hence we investigate the time average 
behavior of the tagged particle's location: 



<52(A) = 



f-A 



t-A 



dt'[xT{t' + A)~XT{tT, (69) 



Bf{n,s) = 



1 



2^/k^ j_ 



dme-v^l"-'"l(/(m,s)) (64) 



since the external force acting on the tagged particle only, 
/(m, s) = 5{m)f{s) (notice that the external force gen- 
erally depends on the time) we find: 



i?/(n = 0,s) 



1 



2y<i 



f{s)- 



for the unbiased case, Fq — Q. For some anomalous pro- 
cesses (5^ ^ (a:^) (e.g. [331). If the time average 5'^ 
is equal to the ensemble average (xip) in infinite mea- 
surement times, the system is called ergodic in the MSD 



sense. First we find ( (52(A) \. There is a relation between 
the correlation function we obtained in the previous sec- 



(65) tions and ((52(A) 



(52(A) 



t-A 



t-A 



dt' [(4(i' + A)) + (4(t')) - 2 {xT{t' + A)xT{t'))] . 



(70) 



Using Eqs. © and ^ for {xrit + A)xT{t)) and (a;|) 
gives: 



<52(A) 



lat 




^3/2 



(t-A) 



1/2 



(2t - A)3/2 3 
2i/2(t_A) ^2172 



Ai/2 + (2 



1/2 



^3/2 



(71) 



for lattice initial condition. When A << t 
= 2a\ 



(52(A) 



lat 



TT 3 f 



When the system is in equilibrium state 



52(A) 



Hence when A << t we find 




(72) 



(73) 



(52) = (52 



where 



p^ = a is the average spacing between particles. We see 
that the ensemble averaged time average MSD is not sen- 
sitive to the way the system was prepared, when A << t. 

On the other hand we find ( S^) ^ (^t) 

/ lat 



lat ' 



A. Relation to fractional Brownian motion 

In equilibrium initial condition, it was shown that the 
problem can be mapped into fractional Brownian mo- 
tion (fBm) with Hurst parameter i/ = 1/4 [H [sUsOJ. 

Hence the process is ergodic in the MSD sense, i.e (<52\ 

is equal to the ensemble average {x'^) and the variance 
of (52 decays to zero when the measurement time is long 
[28j . Therefore, in an experiment with equilibrium initial 
condition, the time average 5^ equals ensemble average. 
Note that the universality of fBm is restricted, since if 
the initial condition are not equilibrium state we get a 
different behavior for the MSD. 



VII. SUMMARY AND CONCLUSIONS 

The generalized diffusion coefficient D1/2 defined in 
(xrity) ~ 2Di/2t^^^ is sensitive to the initial prepara- 
tion of the system. This is a fascinating observation since 
we usually assume that in the long time limit initial con- 
ditions do not influence the asymptotic behavior, and 
that diffusivity is not sensitive to the method of prepa- 
ration of the bath. Moreover we find a sensitivity to the 
initial conditions for two-time correlation function. On 
the other hand in the presence on weak external force we 
find that the two initial conditions give the same result 
for the mean-displacement of the tagged particle since 
it depends on the mean initial conditions. We can con- 
clude that the long-time behavior of the tagged particle 
is determined by the initial condition, in the absence of a 
force field, in the sense that the system always "remem- 
bers" its initial state. Similarly, the time averaged MSD 




FIG. 9: The time average (S'^) when particles are in equilib- 
rium state (upper curve) is compared with the time average 
when particles are initially situated with equal distances with- 
out randomness (lower curve). The simulations sampling is 
representing by the o (lattice case) and * (uniform case), and 
the theory by solid lines (Eqs. (|7H) and (I73|l '). The time mea- 
surement is i = 1000. When A/i << 1 the two ensembles 
give identical results, namely there is no sensitivity to the 
initial condition. 



when A << t is not sensitive to the initial preparation. 
It might be worth while investigating the fluctuations of 
(52 and xt in presence of a field, for the two types of ini- 
tial conditions. A comparison between thermal and non- 
thermal initial condition was discussed recently [40, |41| . 
The elastic model in [40| and the mathematical details 
are different then those considers here but the main con- 
clusion, that initial conditions have ever lasting effect, 
is shared. It seems that non-equilibrium quantities of 
many-body interacting system, like -D1/2 may depends 
on initial conditions, at least in some cases, so it would 
be nice to find more examples to this surprising effects, 
e.g. the dynamics of a tagged particle in one dimension 
evolving according to Hamilton's laws without stochastic 
assumptions beyond the random initial condition [42| . 
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Appendix A: Simulations Details 

As mentioned in the text, in one dimension, a hard 
core cohision event is equivalent to two identical particles 
that pass through each other, and after the particles cross 
each other, the labels of the two particles are switched. 
Instead of relabeling the particles after every collision, we 
let particles pass through each other, and then at time t, 
we relabel our particles. So in fact in the interval (0, t) we 
view the particles as non-interacting. This fact simplifies 
the simulation process by ignoring the collisions during 
(0,i), and sorting the particles labels at time t, and not 
after each collision. 

We choose the system to be continuum in space. Each 
sampling time we move each particle with its probability 
density function: 



P{Xn{t + A)\Xn{t)) 



1 



V-iTTDA 



(x„(t + A)-x„(t))^ 
g 4DA 



(Al) 



where n is the particle's label {n E {—N...N}), and 
Xn{t) is the n-th particle's location at the previous step. 
For each sampling time we sort the particles and find the 
tagged one. For simulations we took 10^ identical sys- 
tems. Each system contains 10, 001 identical particles. 
Notice that the system is finite, therefore we need to con- 
sider times that will be short enough so that the particle 
will not reach the edges, VDt « L, and long enough 
that particles interact with each other a « y/Dt. For 
all the simulations the distance between particles (or av- 
erage distance for the uniform distribution initial state) 
was taken to be L/N — 1. The free particle diffusion coef- 
ficient was D — 1/2. For simulations we used MATLAB 
programming, and used its standard sort and random 
numbers generator functions. 

For the biased case, the tagged particle (the central 
particle) has probability density function 



PixT{t + A)\xTit)) 



1 (x^U + A)-.,U)-i^)^ 



y/AnDA 



(A2) 

instead of the PDF at Eq. (|A1|) . The thermal energy is 
taken to be /csT = 1. 



Appendix B: Normal Distribution of a 

To prove that a = X]i=i ^'^i i^ normally distributed 
when N — > c» we define the kurtosis: 



72 



ik2y 



A^4 - 3^1 



(Bl) 



where the fik is the fc-th central moment of a and is 
defined as: 



^J■k = ((a- {a)f) 



For normal distribution we expect that 72 — > 0. First 
we find the 4th central moment 114. Using the multinom 
formula and the fact that the first moment of San is zero, 

i.e {San - {San}) = gives: 



N 



M4 = 5]((<5«„-(<5«„))') (B3) 

n=l 

+6 Y, (('^a" - {San)f) {{Sam - {6a^)f) . 
For the 2nd central moment we find: 

N 

71=1 

therefore the numerator in Eq. (jBip is: 

N N 

/C4 = X (^{San - {San))j - 3 X [{5a^n) - {San] 



(B4) 



Eq. 



(B5) 



gives: 



1 PRLian)PLL{-an) 
fo„ = { PRL{an)PLR{-an) + PjiR{an)PLL{-an) 

-1 PRR{an)PLR{-an) 

(B6) 
Therefore: 

{San) = (Sal) = PRLina) - PLR{~na) (B7) 

{San) = (Saf-^) ^ PRLina) PLL{-na) + PRR{na)PLR{-na) 

To evaluate k^ we start with : 



N ^ ^ 

E PRLina) = J2 ^^•^(^^) « ^ (B8) 

n^l n— 1 



(generally this term depends on the time t, we take N — >■ 
cx) and t finite). 

Similarly we can prove that ^ Plr oc N,J2 Plr °^ ^ 
etc. hence k^ ex TV. we use the same derivation for the 
denominator and we get A:| c>c N"^ . Finally, we get : 



N AT-i-oo „ 

72 - ^ ^ 0. 



Therefore a is normally distributed: 



„, , 1 (a-(a>) 

P(a) = ^ p 



v/2^ 



(B9) 



(BIO) 



(B2) where Ua = {a^) — 



11 



[1] T. E. Harris, J. Appl. Probab. 2, 323 (1965). 

[2] D. W. Jepsen, J. Math. Phys. 6, 405 (1965). [25] 

[3] D. G. Levitt, Phys. Rev. A 8, 3050 (1973). 

[4] P. M. Richards, Phys. Rev. B 16, 4 (1977). [26 

[5] S. Alexander and P. Pincus, Phys. Rev. B 18, 4 (1978). [27 

[6] S. F. Burlatsky, G. Oshanin, M. Moreau and W. P. Rein- 

hardt, Phys. Rev. E 54, 3165 (1996). [28 

[7] B. Derrida, Phys. Reports 301, 65-83 (1998). [29 

[8] S. Pal et al. J. Chem. Phys. 116, 5941 (2002). 
[9] A. Zilman, J. Pearson and G. Bel, Pyhs. Rev. Lett. 103, [30 

128103 (2009). [31 

[10] E. Ben-Nairn, Phys. Rev. E 82, 061103 (2010). 
[11] O. Flomenbom, Phys. Rev. E. 82, 031126 (2010). 
[12] D. Lucena et al. Phys. Rev. E 85, 031147 (2012). 
[13] D. M. Ackerman, J. Wang and J. W. Evans, Phys. Rev. [32 

Lett 108, 228301 (2012). 
[14] C. Mondal and S. Sengupta, Phys. Rev. E 85, 020402(R) [33 

(2012). 
[15] A. Ryabov and P. Chvosta, J. Chem. Phys. 136, 064114 [34 

(2012). 
[16] P. C. Bressloff and J. M. Newby, Rev. Mod. Phys. 85, [35 

135196 (2013). 
[17] A. L. Hodgkin and R. D. Keynes, J. Physiol. 128, 61-88 

(1955). 
[18] R. I. Macay and R. M. Oliver, Biophysical Journal 7, [36 

545-554 (1967). 
[19] K. Hahn, J. Karger and V. Kukla, Phys. Rev. Lett 76, [37 

2762 (1996). 
[20] Q. H. Wei, C. Bechinger and P. Leiderer, Science 287, [38 

625 (2000). [39 

[21] C. Lutz, M. CoUmann, P. Leiderer and C. Bechinger J. 

Phys. Condens. Matt. 16, S4075 (2004). [40 

[22] L. Lizana and T. Ambjornesson, Phys. Rev. Lett 100, 

200601 (2008), [41 

L. Lizana and T. Ambjornesson, Phys. Rev. E80, 051103 
(2009). [42 

[23] J. B. Delfau et al. Phys. Rev. E 85, 041137 (2012). 
[24] A. Taloni and F. Marchesoni, Phys. Rev. Lett 96, 020601 



(2006). 

L. Lizana, T. Ambjornesson, T. Taloni, E. Barkai and M. 

A. Lomhoh Phys. Rev. E 81, 051118 (2010). 

E. Barkai and R. Silbey, Phys. Rev. E 81, 041129 (2010). 

E. Barkai and R. Silbey, Phys. Rev. Lett 102, 050602 

(2009). 

W. Deng and E. Barkai, Phys. Rev. E7Q, 011112 (2009). 

P. M. Centres and S. Bustingorry, Phys. Rev. E 81, 

061101 (2010). 

P. E. Rouse, J.Chem.Phys. 21, 1272 (1953). 

G. Oshanin et al., in: Instabilities and Non-Equilibrium 

Structures IX, ed. O Descalzi, J Martinez and S 

Rica, (Kluwer Academic Pub., Dordrecht, 2004), p. 33; 

cond-mat/0209611 

E. Barkai, Y. Garini and R. Metzler, Phys. Today 65(8), 

29 (2012). 

Y. He, S. Burov, R. Metzler and E. Barkai, Phys. Rev. 

Lett 101, 058101 (2008). 

B. B. Mandelbrot and J. W. Van Ness, SIAM Rev. 10, 
422 (1968). 

M. S. Taqqu, fractional Brownian motion and long-range 

dependence, in P. Doukhan, G. Oppenheim and M. S. 

Taqqu long-range dependence: theory and applications. 

BirkhSuser (2003). 

A. Taloni and M. A. Lomholt, Phys. Rev. E 78, 051116 

(2008). 

A. Taloni, A. Chechkin and J. Klafter, Phys. Rev. Lett 

104, 160602 (2010). 

C. H. Eab and S. C. Lim, Physica A 389, 2510 (2010). 
L. P. Sanders and T. Ambjornsson, J. Chem. Phys. 136, 
175103 (2012). 

A. Taloni, A. Chechkin and J. Klafter EPL 97 30001 

(2012). 

L. Lizana, M. A. Lomholt and T. Ambjornsson, 

arXiv:cond-mat/1304.1635vl (2013). 

A. Roy, O. Narayan, A. Dhar and S. Sabhapandit, 

J.Stat.Phys. 150, 851 (2013). 



